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^\i ' Abstract 

^ , Nearly-elastic model systems with one or two degrees of freedom are considered: 

the system is undergoing a small loss of energy in each collision with the " wall" . We 
show that instabilities in this purely deterministic system lead to stochasticity of its 

^vq ' long-time behavior. Various ways to give a rigorous meaning to the last statement 

are considered. All of them, if applicable, lead to the same stochasticity which 

^> , is described explicitly. So that the stochasticity of the long-time behavior is an 

P I ' intrinsic property of the deterministic systems. 
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T^ ■ 1 Introduction 

C^ , Consider the one-dimensional motion of a unit-mass particle in a smooth potential 

fv^ . F'iQ) in an interval [01,02] with elastic reflection at oi and 02. Let F'[ai) > and 

F'(o2) < (Fig. 1(a), the potential F{q) is shown there as the bold-face line): if g G 

[01,02) and p > or g G (01,02] and p < 0, the trajectory starting at (g,p) moves 

/\ . according to equation q{t) = —F'{q{t)); the trajectory jumps instantaneously from {ai,p) 

Cd . to (oj, —p), if i = 2 and p > or if i = 1 and p < 0. If the initial velocity is large enough, 

the particle hits both "walls" oi and 02 and performs periodic oscillations shown in 

Fig. 1(b). Let now the walls be not absolutely elastic. It is natural to assume in certain 

situations that the loss of energy is proportional to the speed at the collision point 

with the wall: if the particle hits Oj with a speed v, it is instantly reflected with the 

speed —v{l — eci{v)), i £ {1,2}. Here Ci{v) are positive smooth functions, < e << 1. 

The coefficient (1 — eci{v)) is called the coefficient of restitution. Denote by g'^(t) the 

position at time t of the particle performing nearly-elastic motion, p^{t) = q^{t). It is 
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Fig. 1: The 1-dimensional mechanical model 
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clear that, for each e > and t large enough, g'^(t) will be situated in an arbitrary small 
neighborhood of one of the points ai,a2, 03 (We assume that the potential F(q) has the 
form shown in Fig. 1(a) and 03 is the unique local minimum inside [01,02]). 

To be specific, assume that the initial speed q^{0) = p^(0) = vq is large enough 
: Vq > Y^2maxF((7). Then q^{t), for e small enough, hits both points ai and 02, but 
lim g^(t) = q^ exists and q^ = ai or q^ = 02. Which of these two points appears as 
lim 9^(t) depends on e in a very sensitive way so that limg^ 

t-s>oo e4,0 

exist. We show that the final position q^ converges, in a certain sense, as e 4, to a 
random variable distributed between the points oi and 02. There are various ways to 
make the last statement rigorous. But they all lead to the same distribution between ai 
and 02, so that the stochasticity of the system as e << 1 is an intrinsic property of the 
system. 

The perturbed system q'^it) for < e << 1 has fast and slow components. The 
fast component consists of the motion along the non-perturbed trajectory. To describe 
the slow component, consider the graph F obtained after identification of points of 

each connected component of every level set of the Hamiltonian H{q,p) = \- F{q) 

(Fig. 1(c)). Denote by y : n — )■ F the identification map of the phase space n = 
{(?) p) £ I^^ • oi < (7 < 02} of our system on F. The slow component of the motion is 
y((7^(t), q^{t)) (compare with [6], Ch.8). Number the edges of the graph (F = /i U /2 U 
... U I5 in the Fig.l(c)). Then Y{q,p) = {H(q,p), K{q,p)), {q,p) G Fl, where K{q,p) is 



p2 

the number of the edge containing Y{q,p) and H{q,p) = h F{q). The pair {H,K) 

form a global coordinate system on F, and the slow component is the motion Yf = 
{H{q%t),q'{t)),K{q%t),q%t))) on T. 

We prove that, in a certain sense, the rescaled slow motion Y^, converges weakly 
to a stochastic process Yf on T. Inside the edges, 1^ is a deterministic motion and the 
convergence follows from more or less standard averaging principle (see, for instance, 
[1], Ch.lO). The stochasticity appears due to a branching at the interior vertices of T 
(compare [2]). 

Note that in the case of system shown in Fig.l, the phase trajectory {q^ (t) , q'^ {t)) 
with (Jq = vq > /2 max F{q) never enters the domain {{q,p) G n : Y{q,p) G I^}, so 

V ai<q<a2 

that if the damping of the system occurs just on the walls oi and 02, the stochasticity 
of the limiting slow motion is concentrated at the vertex corresponding to the absolute 
maximum of the potential. Therefore we will consider in more detail the case when F{q) 
has just one local maximum on [01,02]. If the potential has many wells and the system 
is losing energy not just at the walls ai, 02, one should take into account the whole 
graph r even if the initial speed is large. 

In the next two sections we consider a model problem where the limiting slow 
motion inside the edges is just a motion with constant speed. This system has some 
damping not just in the ends of the interval [01,02] and models systems with multi- 
well Hamiltonian. To give a rigorous meaning to the statement that the slow motion 
converges to a stochastic process, we, first, perturb the system g^(t) by a stochastic 
perturbation of a small intensity 5 << 1. Then in this already stochastic system q^' (t), 
we consider the slow component Y^' = Y{q^' {t),q^' (t)). We show that Y^,' converges, 
as first e J, and then (5 J, 0, in a certain sense to a stochastic process Yj on T. This 

double limit limlimK^' = Yt exists and is the same for a broad class of stochastic 
510 £4-0 

regularizations. It is in this sense that the slow component of the purely deterministic 
system q^{t) should be approximated by the stochastic process Yt on T. 

In section 4 we come back to the problem mentioned in the beginning of this section. 

Of course, one can consider similar questions for systems with more than one degree 
of freedom. Let, say, F{q), q £ M? he as shown in Fig. 2. The particle of a unit mass 
moves inside a convex domain G with a smooth boundary dG according to the equation 
qt = — VF(gi) and undergoes an instantaneous mirror reflection on the boundary. If 
the collisions with the wall dG are absolutely elastic, the particle will move forever on 

the energy surface H{p,q) = — — \- F{q) = Hq. But if the particle loses energy in 
each collision similar to the case of one degree of freedom, it can, eventually, be found 
near one of the local minima of F{q) on dG (points Oi and O2 in Fig. 2). We consider 
many-degrees-of-freedom problems for model potentials in Section 5. 




Fig. 2: A 2-dimensional problem 

2 One-degree-of-freedom Model Problems 

Consider a model of the system with several potential wells. Suppose a parti- 
cle of unit mass moves freely in an interval [qi , g„] with elastic reflection at the ends 
of the interval if the initial velocity is large enough. Let a finite number of points 
q2,q3, ■■■, Qn-i G (<?i) <?n) and < p2 < ••• < Pn-i are given such that if the particle starts 
with a velocity po at a point qo G [qi,qn] then it moves freely in fe_(po,go)'9i+(po,go)] 
with instantaneous reflection in the ends of this interval, where i_ = i-{pQ,qQ) and 
i+ = i+{po, qo) are defined by the conditions (see Fig.3): 



Qi- < qo< qi+ , Po < mm{pi_,pi_^_) , po> pk for k £ {i : q^ £ iqi_,qi+)}. 

T) 

The energy — is preserved in this system. This is an approximation of the motion in 
a potential which has sharp maxima at points qi and is close to a constant between qi 
and qi+i- Just for brevity, we assume that the walls between the maxima have the same 
depth. 

Assume again now that the collisions with the walls are not absolutely elastic: if 
the particle hits the wall at the point g^ with a speed p, \p\ < Pk, it is reflected from qi^ 
with the speed —p{i — eckip)), where Cfc(-) are defined as before but now (1 — eck{p)) 
represents the coefficient of restitution for the wall {{q,p) G M^ : < p < Pk,Q = Qk}, 
and < e << 1. Then for each e > 0, the velocity tends to zero as t — )■ oo. On each 
velocity level Pk, Pk < POi phase trajectory goes to one of the two "wells" separated by 
qk- Corresponding phase space Fl = {{q,p) £ K^,9i < q < q-n} is shown in Fig.3 for 
n = 4. Notice that actually each interior wall in the phase space n consists of two sides 
(the broken line represents another side in Fig.3). For simplicity we will identify the two 
sides in our investigation and we assume that the coefficients Ck{-) is the same for both 
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Fig. 3: l-diniensional model problem 



faces of the wall q = Qk- 

We will see that the "choice of the well" at each q^ is very sensitive to e as e J, 0, 
and the behavior of the trajectory becomes, actually, stochastic in a certain sense. Now 
we will give the exact meaning to this statement and describe the limiting stochastic 
process. 

Let us, first, consider the case when the model system has only two trapping wells. 
An example of such a model is shown in Fig. 4. A particle with unit mass starts its motion 
from a point x = {qo,Po) in the phase space n and has instantaneous reflection each time 
when it hits the boundary q = —ai (corresponding to well £i) ov q = 02 (corresponding 
to well £2)- Each time when the particle hits the boundary q = —ai or g = 02 it 
will be reflected and move to the point (— 1,— j)(l — eci{p))) or (1, — p(l — ec2{p))), 
respectively. The coefficient of restitution for the middle wall (notice that it has two 
sides) (7 = is 1 — ecs^-). After the particle enters one of the wells £1 or £2, it hits 
corresponding boundaries and loses energy according to the coefficient of restitution of 
the corresponding boundary. We denote such a motion by Xp^ = Xf = {q^{t),p^{t)), 
where (9'^(0),p^(0)) = ((^OiPo) = a;. The superscript x represents the starting point 
of the motion. Here and henceforth, when we use the notation Xf, it means that 



we have omitted the superscript x. Put H^{t) 
X 



We rescale time and let 



^ — X^i. As before, we consider the slow component of Xf which is the projection 
{H{XI),K{XI)) of Xf onto the graph F. Here H{Xf) = H%t/e). Since the 
coefficient of restitution depends only on the magnitude of velocity but not its direction, 
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Fig. 4: The case when the model system has only two wells 



i.e. Ci{—p) = Ci{p), i = 1,2, we can write Ci{—p) = Ci{p) = Ci{\j2H). Now we prove: 



Lemma 2.1 Within each edge of the graph T, as e \. 0, the first component of the 

process Y^ , i.e. H{Xf) = , converges uniformly to a deterministic motion 

H{t) which satisfies the differential equation 

dH ci{V2H) + C2iV2H) 



dt 



n(H) 



H , H{0) = Ho on I3 



and 



dH ClUm+CsiVmrr J 

-df = -' T^iH) ' ^' 

d_H ^ c,iV2H)+c,iV2H) 
dt T2{H) 



(2.1) 

(2.2) 

(2.3) 

are the corre- 



respectively. Here n{H) = ^^^L±^^ Ti{H) = ^, T^iH) = -^ 

Vm Vrn V2H 

sponding periods of motion on phase picture n along the energy level H for each edge of 

r. 



Proof: Let us prove, for example, (2.1). The proofs of (2.2) and (2.3) are exactly 
similar to the proof of (2.1). We prove this result by using a slight modification of the 
standard method of justification of the averaging principle (compare with [1], Ch.lO). 
First of all, we view the whole contour (see Fig. 4) A ^ B ^ B' ^ A' ^ A on which 
the particle has a fixed amount of energy H as a loop Gh (i-e. identify A with A' and 



6 



B with B'). One can introduce coordinates (Q^, ,H^, ) to describe the motion of the 
particle on this circle. Here Q^/^ = qf/^ + 2(fc — l)(ai + 02) if the particle is the fc-th 
time on the upper half plane of the phase space, and Q^,^ = —qfi^ + 2(fc — l)ai + 2ka2 
if the particle is the k-th time on the lower half plane of the phase space. We identify 
Qw with Qw +2k{ai +02), k gZ. The motion (Qw ,H^, ) has a fast component Qw 
and a slow component H^, and they satisfy the following equation 



Qt/e = -J^H^/s 



H 



tie -= -iMQt/e - «2) + A(Qf/^ + ai)](2c(Qf/,, ^2H^/^) - ec^Qt^^, ^'^Hl^Ml/, ■ 

(2.4) 



Here A(-) is the Dirac (5-function and c(Qf , , .r2H^jJ = ci( /2iJ|, ) if Q 



tje' 



tje' 



'tie 



-a\ 



and c(Qf/^, ^2Hli^) = c,{^2Hli^) if Qf/^ = a,. 

Now we introduce an auxiliary function u = u{Q, H) satisfying the following dif- 
ferential equation 



■du 



'2H—- = A-A, QgGh 
oQ 



(2.5) 



where 



A = -[A(Q - 02) + A(Q + ai)](2c(g, Vm) - ec\Q, V2H))H, 



and 



A 



^ci{V2H) + C2{V2H) - e{cl{V2H) + cl(V2ff)) 

mm 



is the average of A along the loop A^B^B'^A'^A. 

One can check the following properties of u{Q,H): (i) the differential equation 
(2.5) describing u really has a well-defined solution u{Q,H) on the loop Gh', (ii) the 
solution u{Q,H) is uniformly bounded together with its first derivatives for Q G Gh, 
H £ [H{0),M] for any M G {H{0), 00). 

Now, by the fundamental theorem of calculus, we have: 



HQtle^H^le)-^iQo,m) 

du ■ f du ■ 

- I {A-A)ds + 0{1). 
£ Jo 



Therefore we see that max 

0<t<T 

Now, we compare the two differential equations: 



{A - A)ds 



< Ce for some positive C > 0. 



% = -[A(Qf/,-a2)+A(Q^/,+ai)](2c(Qf/,, ^2Hl^^)-ec\Ql/^, ^2i^))if,% , HI, = Ho; 
Let m(T) = max {Hf, — H(t)\. We have, that 

0<t<T */^ 

\Ht/,-H{t)\ 

I {AHl, - AH{s))ds + 2eB 
Jo 

< [ \A- A\Hl,ds + / A\HI, - H{s)\ds + 2eB. 
Jo Jo 



The last inequality implies that m{T) < [CD + 2B)e + E m{s)ds, where B is 

Jo 
„2 (' ^/2ff) -\- c^l \/2H) 

the uniform bound for ; — , D is the uniform bound for H"^ , , and 

_T3{H) '/^ 

E is the uniform bound for A. By using Gronwall's inequality, we see that m(T) < 

(CD + 2B)ee^^, i.e. lim max \Hf, - H(t)\ = 0. ■ 
£4-0 0<t<T' */^ 



Notice that in this case Ti{H) ~ H^2 as -ff — )■ H{0), this means that Yf will enter 
the interior vertex O of the graph F in finite time. Thus one might ask what is the 
behavior of the motion Yf at the interior vertex O. First, we regularize this problem by 
considering small stochastic perturbation of the initial conditions in our model. 

Lemma 2.2 Let U{x,6) = {y £ n : \x - y\ < 6} , x £ n, 6 > 0, Uf = {x £ n : 
X^'^ eventually enters the well Ei\, i G {1,2} (see Fig. 5), where X':'^ is the perturbed 
trajectory starting at Xq'^ = x £ H. Assume that Hq = H{x) > H{0). Then 

n{Ul^U{x,5)) ^ ci{H{0)) 
Sio eio fi{UI f]Uix, 6)) C2{H{0)) ' 

where pi is the Lebesgue measure in M?. 

Proof: Without loss of generality one can assume that pq > 0. We see, from Fig. 5, 
that U{x, 5) is covered by narrow shaded or white strips, where shaded strips belong to 
Uf, and white strips belong to Z^l- For a fixed x, consider the nearest to x shaded and 
white strip (ordered as a pair, shaded on top). The upper p-level of the shaded strip 
of the pair is denoted by a. And the lower p-level of the shaded strip (also the upper 
p-level of the white strip in the pair) by b. Let the lower p-level of the white strip in the 
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Fig. 5: Regularization by stochastic perturbation of the initial condition 

pair be c. Let n = n{e, x) be the number of collisions that the particle made with either 
of the walls q = —ai or q = 02 before entering one of the wells. 

Let us first estimate the p-width of the shaded and white strips. Put 

A{a) = a{l - eci(a))(l - ec2(a(l - eci(a)))) = a(l - e/f (a) + eVl(a)) 

where /f (a) = ci(a)+C2(l— eci(a)) and /2(a) = ci(a)c2(l— eci(a)). Then the boundaries 
of the next shaded and white strip has p-height A{a) = c,A{b),A{c), respectively. 
We see that 

\c-a\ = \A{a) -a\= ea|/f (a) - e/|(a)|, 

So that 

Mie <\c-a\< M2£ 

for some Mi,M2 > 0. Therefore, n = n(x,e) ~ 0{—). 
We have 



A{a) - A{b) 

= a-b- e{afl{a) - bfHb)) + £'(«/!(«) " bfl{b)) 

= a-b- £{ul{b){a - b) + u|(0(« - bf) + e^K(6)(a - b) + vl[ri)[a - b) 

= {a - b){l - eulib) + 0{e^)) 

where Ui, Vi, i £ {1, 2} are bounded smooth functions, and C, 77 G [b, a]. 
Similarly, 

A{b) - A{c) = (6 - c)(l - eulic) + 0{£^)). 



These inequalities imply that 
A{a) - A{h) 



A(h) - A[c) 

a-b { I - eu\{h) + 0{e^] 



b-c Vl-enf(c) + 0(e2 

a-b f e{ul{b) - u\{c)) + 0{e^) 

b-c V l-enf(c) + 0(e2) 

{l + 0{e^)). 



?LzA(^^n(r^ 



Thus after n = n{e,x) r^ 0{—) steps the change of ratio of the p-width of the 
nearby shaded and white strips is asymptotically (1 + ©(e^))"^^'^) ~ 1, and our lemma 
then immediately follows from the definition of the sets U{x,6), Uf and Wf- ' 



In order to include our result in the framework of weak convergence of stochastic 
processes in the space of continuous functions, we still need the following technical 
construction. Notice that, for fixed e > 0, the function H^{t/e) as a function of t, is a 
step function. We can assume that it is continuous from the right and having limit from 
the left. We now connect neighboring discontinuity points of the graph of H^{t/e) by 
straight line segments. The resulting trajectory is denoted by H^{t). For fixed T > 0, 
we see that H^{t) G Cqt- For a positive constant C > and all < t < T, we have 

\H%t/e)-H%t)\<Ce. (2.6) 

Since the time distance between two discontinuity points is bounded from below 
by the inverse of the initial velocity, the slope of all the line segments of the graph of 
H^{t) is bounded uniformly in e G (0, 1]. Thus the family of functions {i:f^(t)}e>o is an 
equicontinuous family in Cqt- Also, it is uniformly bounded for < t < T and e G (0, 1]. 
Thus by Ascoli-Arzela Theorem the family {i7^(t)}£>o is compact in Cqt with uniform 
topology. 

Let ^s be a two-dimensional random variable with a continuous density fs{x) such 
that fs{x) > for |a:| < 5 and fs{x) = for |x| > 6. Let X^' = (pf ,ql' ) be the 
trajectory of the nearly-elastic motion starting at Xq = x + ^s (the point x = {qo,Po) G 
n is fixed, and we excluded it from the notation). Put X^' = X^', . 

Consider a stochastic process Yt = {H{t), K(t)), < t < T, on the graph T defined 
by the conditions: Yq = Y{qQ,pQ); H{t) is a deterministic motion inside each edge of 
r satisfying equations (2.1)-(2.3) respectively; if Yq = ^(qojPo) £ ^3) the trajectory Yt 
reaches the interior vertex O G F in a finite time, instantaneously leaves O and enters the 
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1 1 .1 ci(H(0)) C2{HiO)) 

edges i 1 or I2 with probabilities pi = — , pi — 



ci{H{0)) + C2{H{0)) ' "^ ci(//(0)) + C2{H{0)) ' 
respectively. The conditions listed above define the process Yt in a unique way (in the 

sense of distributions) . 

Taking into account (2.6), the remark concerning compactness and Lemmas 2.1 

and 2.2, we get the following 

Theorem 2.1 The process Y^' = {H(X^' ),K{X^' )) on T converges weakly in 
the space of continuous functions [0, T] — t- F equipped with the uniform topology to the 
process Yt = iH{t), K[t)) when, first e 4, and then 5 4, 0. 



We now discuss the case when the model problem has more than two energy trap- 
ping wells. It turns out that in the multi-well case, the perturbation of the initial 
condition, in general, will not lead to a regularization of the problem: the limit of Y^' 
as e 4, may not exist. 

Consider an example as shown in Fig. 6. There are 3 trapping wells, 1, 2 and 3. 
Wells 2 and 3 are separated on a p-level p^ and they combined together are separated 
from well 1 on a p-level p*. Pa > pI- We call the wells 2 and 3 combined together well 
4. Well 1 and well 4 combined together is called well 5. Suppose the final strip A as 
drawn in the picture is for those initial values that finally enter well 4. The boundary 
restitution coefficients are positive constants ci and 02 = 03 = 04 = 05 = 0(7^ ci) as 
shown in Fig. 6. 

For each e > 0, the shadowed strip A on the top of wells 1 and 4 is moved by the 
nearly-elastic dynamics to the well 4. The strip A has the width ce. The well 4 between 
p-levels pI and p* is also covered by strips of width ce (suppose that e is such that an 
integer number of such strips is situated between p^ and p*). The highest strip B in the 
wall 4 is moved by our dynamics (as a whole) alternatively either to well 2 or to well 3 
as e 4, 0. But the strip A in one step (one reflection from the right exterior wall) is going 

to the strip B. Therefore, for the sequence of e 4- for which — is an integer, the 

e 
nearly-elastic trajectory starting from any x = {qo,Po), Po > P*ai fo^ ^ large enough is 

distributed alternatively either between walls 1 and 2 or between walls 1 and 3. This 

means that the limiting distribution for Xl' as e 4- does not exist, and the problem 

cannot be regularized by random perturbations of the initial point. 
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Fig. 6: The case when the system has more than two trapping wells 

3 Regular izat ion by stochastic perturbation of the dynam- 
ics 

Let us now consider small stochastic perturbation of the dynamics rather than the 
initial condition. We will show that, this regularization works when the system has any 
number of trapping wells. But let us first start from the case of system with only two 
trapping wells. 

The particle of unit mass starts its motion from x = {qo,Po) with po > (see Fig.7). 
Each time when the particle hits the wall q = —ai or q = a2, it instantaneously moves 
to (— ai, — p(l — eci — e6r]f^)), p < and (02, — p(l — £C2 — ed^k)), P > 0, respectively, 
and then it is reflected. We assume for simplicity that ci,C2 are positive constants, 
< e << 1 and < 6 « 1. The restitution coefficient for the interior wall is assumed 
to be 1 — ecs — e5C,k for both faces. And C3 is also a positive constant. 

Here ^fc, %, C/o ^ = 1)2,..., are independent sequences of i.i.d random variables 
with continuous densities; the random variable with subscript k is used in the restitution 
coefficient when the trajectory hits the corresponding wall fc-th time. In the following 
we also use the general notation ^, 77, C, to denote independent r.v.'s which has the same 
distribution as ^fc, rj^ and Cfcj respectively. We assume for brevity that random variables 
^, r] and (" are bounded with probability 1. Then, without loss of generality, one can 
consider just strictly positive ^, rj and Q: P{a < ^ < /?} = P{a < r] < 13} = P{a < 
C < /3} = 1 for some < a < (3 < 00. Actually, one can replace the boundedness of 
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Fig. 7: Regularization by stochastic perturbation of the dynamics 

random variables by their positivity and the finiteness of their moments (together with 
the existence of a continuous density) . 

The resulting process is denoted by X^' = {q^' {t),p'^' (t)). We rescale time and 



put Xf = Xlf^. Let Yt' 



{H{X^' ),K{Xp )) be the projection of X^' onto the 



graph r (see Fig.7). Put H^'^{t) = H{Xl' ). For fixed e > 0, this is a stochastic 
process with jumps and each sample path is a step function which is right continuous 
and has limit from the left. Using the same construction as in the end of Section 
2, we obtain from H^' {t) a continuous piecewise linear process H^' {t). We also put 



>£,<5 



-£,<5\ 



?e,5 



Y^ ' = {H{X^' ),K{X^' )). Since ^ and rj are bounded, 

P{\H''^{t)-H''^{t)\ <Ce) = l 

for some constant C > 0. 

Now we are going to prove the weak convergence of Y^ ' in the space CqtO^) of 
continuous functions [0, T] — )■ F provided with uniform topology, as first e J, and then 
(5 4, 0, to a stochastic process Yt on the graph F. We do this through a series of lemmas. 



Lemma 3.1 For fixed 6 > 0, the family {Y^' {t)}^yQ is tight in the space Cot(F). 



Proof: We have to verify the following (see [7], Ch.6, Theorem 6.4.2): there exist 
a, /3 > 0, s.t. for any /i > 0, any 0<t<t + h<T, any e > 0, 



E p(y/,y/+,) <Mh 



,l+/3 
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Here p(x, y) denotes the distance between the points x,y £ T defined as follows: 
if X and y lies in the same edge Jj (i = 1,2,3) of the graph T, then p{x,y) is just the 
usual Euclidean distance between points x and y within the line segment; otherwise, 
suppose X £ li and y £ Ij, i ^ j, then consider a path lying on the graph F connecting 
X and y, we denote p{x,y) be just the minimal length of such paths. Since our ci, 
C2, C3, and ^, ?7, C are bounded with probability 1, and the period T{H) of the non- 
perturbed (elastic) motion is bounded for < i:f < 00 and separated from zero, we 
see that P{p{Yf,Yf,j^) < Mh) = 1 for some constant M > 0. Therefore one can take 
a = 2,/3 = 1, and the statement follows. ■ 

Lemma 3.2 Let H'^'^{0) = Hq > H{0). Within each edge of the graph T, as 
e 4, 0, the process H'^' (t), for < t < T < 00, converges uniformly in probability to a 
deterministic motion H (t) which is defined by the equations 

H\t)=(''-l±^^l±^^^^±^ + HM'\ onI„ 



and 



t-toCi+C3 + 5{EC + E.T]) ^ ^,^,_i/2 



H\t) = ( ^" ^ al + H{0)-'^' ] , onh, 



-2 



t - to C2 + C3 + S{E^ + EC) 



H\t) = --;0^2^C3^.,^,^^^; ^ F(0)-V2 ^onl, 



-2 



\/2 a2 

respectively. Here H{0) is the energy corresponding to the interior vertex O and tg 



Cl + C2 + SE{^ + T]) 

vertex O. 



is the time for the motion H (t) to come to the interior 



Proof: The proof of this lemma is similar to the proof of Lemma 2.1 and we use 
the same notations. The system (2.4) should be replaced by the following system: 



1 



Ht,e = -[^{Ql/e - «2) + A(Qf/, + ai)](2c(Qf/J + 25 ■ riQ^^ - eb{Qip)HI^^ . 

(3.1) 

Here A(-) is the Dirac (5- function so that the right hand side of the last equation 
in (3.1) is not zero just if Q^, = oi or Q^, = 02; ^(Qwg) = Ck and c{Ql,J = C2 
when c{Q^,J = 02; r{Q^,) = rj^ and c{Q^,J = c\ when c((5wg) = a\. The function 

KQt/e) ~ ( ^iQl/e) + ^ ' ''(^t/e) ) ^^ Uniformly bounded for < t < T with probability 



one. After that, using the same arguments as in Lemma 2.1 and the law of large numbers. 
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we obtain an equation for the limiting slow component H{t) on each edge. For constant 
Cj, these equations can be solved explicitly and we get the statement of Lemma 3.2. ■ 

Lemma 3.2 implies the following: 

Corollary 3.1 Within each edge of the graph T, for the time < t < T < oo, 
the process H'^'^{t), as e,6 I 0, converges in probability to a deterministic motion H{t) 
defined by the equations 

Hit)=(l=2l±£l + HM'\onI, (3.2) 

V V 2 oi + 02 / 

and 

Hit) = (^^ + HiOyy^^ "' , on h (3.3) 

Hit) = (^^ + HiO)-^-) " , on h (3.4) 

respectively. Here to = ^ 



Cl +C2 

Let us consider now the slow motion near the interior vertex O. We first prove the 
following auxiliary lemma concerning random walks. Let {^fc} and {f?fc} be independent 
sequences of i.i.d random variables. Assume that the random variables have continuous 
densities and P\a < ^ < /3} = P\a < ?? < /3} = 1 for some < a < /3 < oo. Put 



Dq — X , S2m — X + 2_J^^^ ^ ^^' ' "^Sm+l — 3^ + '52m + Cm+1 • 
fe=l 

It is clear that S"^ = x + S*^, P{na < S^< n/3} = 1. 

Define Tn as the first time m when S^ is greater than n\: Tn = min{?TT, : Sf^ > 
n\}. 

Since E(,^fc + 77^) > 0, the law of large numbers implies that P{Tn < 00} = 1 for 
any x G M\ A > 0, n G Z. 

We use two equivalent types of notations: The initial point x of the random walk 
will be included either as a superscript, like S^, Tn , or as a subscript in probabilities 
and expected values so that KxfiSn) = IE/(5^), Px{Tn < t} = P{Tn < t}. 



Lemma 3.3 Under the conditions mentioned above, lim -Po{t„ is even} 

n— >oo 



lim Po{t^ is odd} 



,^ + E7?' n^oo "^ " ' E^ + E?7 
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Proof: Put mi(n) = mi 



nX 






if the latter integer is even, and mi{n) = rrii 
nX 



1 otherwise. It is clear that S^_^ < -—. Put M = K{$,k + Vk), D = Var{$,k + Vk)- 



mi 



Let mi(n) = mi = — ^. Let fmi{x) be the density of S^^ — miM. 

It follows from the central limit theorem that for each 5 > one can choose N 
N{6) so that 



P{-N^Jmi{n) < S^^ - miM < N^/^iin)} >l-d 

for all large enough n. 

Using the Markov property of the random walk S"^, we get 



(3.5) 



Po{Tn is even} = KoPs^^^^^{t^ is even}. 
One can conclude from (3.5) and (3.6) that 



(3.6) 



Po{t,^ is even} 



NVm 



NVm 



1 

_/mi(a;)P^+miM{7"^ is even}dx 



<5. 



(3.7) 



Since our random variables are bounded and have density, one can apply the local 
central limit theorem to S**^ / ■,: 



J mi \X) 



1 



\/27rmiD 



exp( 



2miD' 



< 



\fm[ 



(3.8) 



1 {k\ (k) 

uniformly in a; G M , here and later, we denote by p) such sequences that lim p) = 0. 

^— s>oo 



Divide the interval 



into 



\Anii 



n) 



equal intervals Ii 



-Ny^mi{n), Ny^mi{n) 
If. It is clear that r = r{n) ~ n^'^ and the length |/fc| of each interval Ik is of order n^'^ 
as n — 7- oo. Because of (3.7) and (3.8), 



r{n) 



Pq{Tj^ is even} — >^ / exp 



fe=i'''^'= 



x^ -.Px+miMiTn is even}dx 



2rniD 



\/2TTmiD 



<5 + pi2). (3.9) 



On each I^, k G {1, ...,r}, choose a point z/. such that 



exp(- 



2miD 



)dx = |4|exp( 



2m^D 



li X ^ Ik., then 



exp( 



2miD' 



exp( 



2miD' 



< pI^^ exp(- 



2miD' 



(3.10) 
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where lim /0„ = since |x| < N\/mi{n), \zj^\ < N\/mi{n), \x — Zk\ < Const • n^''^, 
mi{n) is of order n as n — )■ oo. We conclude from (3.9) and (3.10): 



<6 + p(^K 
(3.11) 



Po{t^ is even} - Vexp(- J' )— =^= / P^+^,m{t-^ is even}-— 
^ 2miD ^j2TimiDJi^^ |4| 

Note that, since the random walk Sn is invariant with respect to shifts, 

Px+An{Tn is even} = Pa;{r^+^ is even}. 
Using this fact and taking into account (3.8), we conclude that 

/ i-j-^Px+miM{rn is even} = / i^-r P^Jt^' is even} + -^ (3.12) 

for an appropriate constant A'. The integral in the right hand side of (3.12) is nothing 
else but the expected value E of the random variable X uniformly distributed on J^. 
Let Xx be the indicator of the set {tu is even} in the sample space. Then 

/ ^P.{r„^' is even} = Wxf = EExJ^^ . (3.13) 

For given sequences {^fc}, {f]k}-> 

g w ^ Vi+m + - + ri.^ + ^n 

where f„ and k„ are random variables such that lim t'n = oo and |k„| < /5. 
Using the strong law of large numbers, we conclude from (3.14) that 

lim Ey'^ = — , o.e. . 

n^oo ^^ E^ + Er/ ' 



Since |Exj^ | < 1, the last equality and (3.13) imply that 



f dx ^ , y . , 

lim / T-r-rPx\Tn is even} 



E77 



E^ + E?? 
From (3.11), (3.12) and (3.15), we derive: 



(3.15) 



Er; 
Poi'^n ^-5 even} 



<25 + p^^\ 



E^ + Er/ 

This bound implies the first statement of the lemma. The second statement follows 
from the first. ■ 
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Now we have, for our process X^' , 



Lemma 3.4 

Cl 



lim lim PjXf ' finally falls into the well E\ | = 

540 £40 "^ * -^ ^•' ^J C1+C2' 

limlimPjXf' finally falls into the well £0} = ■ 

Sio eio ^ t J y J ^J C1+C2 

(see Fig. 7) 

Proof: When the particle coUides with the wall q = — ai, the absolute value of its 
velocity changes as follows: 

\Pnew\ = |Po/d|(l - eci - e5^k)- 
Similarly, when the particle collides with the wall q = 02'- 

\Pnew\ = |Poid|(l - eC2 - £(5%). 

Take logarithm of these two equalities: 

In \pnew\ = In \poid\ + ln(l - eci - ed^k), 

In \pnew\ = In \poid\ + ln(l - ec2 - e5r]k), 

respectively. Thus we can view the problem of determining whether X^' enters well £i or 
£2 as the following random walk problem. We start from In |po|) and alternatively jump 
backward with steplength — ln(l — eci — e5^k) (at step 2k — 1) and — ln(l — ec2 — e^rjk) 
(at step 2k). If the random walk finally jumps over In |p(0)| at an odd step, then X^' 
enters well £i] otherwise X^' enters well £2- 

Now we further simplify the problem: we start from 0, and alternatively jump 
forward with steplength C/^' = — ln(l — eci — eS^k) (at step 2k — 1) and V^' = 



ln(l — ec2 — £<J?7fc) (at step 2k). If the random walk finally jumps over — In 

an odd step, then Xl' enters well £1] otherwise Xl' enters well £2. 



Po 



p{0) 



at 



Put n 



A = ln 



Po 



. Taking into account that U^' = (ci + S^k) + 0{e) 



p{0) 
and Vf' = (c2 + (5%) + 0(e) as e 4, and applying Lemma 3.3 we get: 



e 5 ■ -1 ■ ■ ^ k 



limlimP{Xj' finally falls into the well £1} = lim lim 



EU'/ 



* ^in pin WTT'^'^ I TI7T/£'<5 



510 £40 " 510 £40 Ef/^'" + ^V^' Cl + C2 
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The second statement of the lemma follows from the first. ■ 

Let Yt = {Ht, Kt), Hq > H{0), Kq = 3, be the process on F defined inside the edges 
by formulas given in Corollary (3.1); when Yt reaches at time to = {H{0)^ ' — H^ ) 

C\ ~r C2 

the interior vertex O, it goes immediately to /i or I2 with probabilities pi = and 

Cl +C2 

P2 = respectively. 

Cl +C2 

Combining Lemma 3.1 and Corollary 3.1 and Lemma 3.4 we have the following 

Theorem 3.1 The process Y^' = {H{Xl' ),K{X^' )) on T converges weakly in 
uniform topology in the space of continuous functions [0, T] — )■ F as first e J, and then 
(5 4, to the stochastic process Yt = iH{t), K{t)) on F. 

Remark: As we have seen in Section 2, in the case of two wells, the problem can 
be regularized by stochastic perturbations of the initial conditions. Theorem 3.1 shows 
that the regularization by perturbations of the dynamics lead to the same limiting slow 
motion 1^ on F. 

Now we consider the case when the model system has more than two trapping wells 
(see Fig. 8). Then the regularization by perturbation of the initial conditions, in general, 
does not work. Let the coefficients of restitution of the left and right boundaries of the i- 
th well be 1— ec^ —eSS,^ for the left boundary and 1— ec2 ~^^'% ^°^ ^^^ right boundary. 
In Fig. 8 i G {1,2,3,4,5}. A random variable with the subscript k is used when the 
particle hits a wall the /c-th time; c]^ , C2 are positive constants; {Cfc IfcLi) {Vk iT=i ^'^^ 
sequences of i.i.d random variables. We write c^*-*', rj^"^' to denote independent random 
variables that have the same distributions as any of the random variables from the 
sequences {S,f}t=i, {^i'^l^i- We assume that P{a < ^(*) < (3} = P{a < r/« < /3} = 1 
for < a < /3 < 00. 

The resulting process is denoted by X^' = (g^' (t),j)'^' (t)). We rescale time and 
put X^' = X^]. Let y/' = {H{Xl' ),K{X^' )) be the projection of X^' onto the 
graph F (see Fig.8). Put H'^'^it) = H{X^'^). For fixed e > 0, this is a stochastic 
process with jumps and each sample path is a step function which is right continuous 
and has limit from the left. Using the same construction as in the end of Section 
2, we obtain from H^'°{t) a continuous piecewise linear process H'^' [t). We also put 
Yt' = iH{xf),Kixf)). 

The averaging procedure is the same as before. One can see that within each edge 
/j, as first e \, and then 5 \,0, the process H'^' (t) converges to a deterministic motion 
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Fig. 8: The case when the system has more than two trapping wells 



H{t) which satisfies the difi'erential equation 



dH 
It 



Cji-C2 



(3.16) 



where i G {1, 2, 3, 4, 5}. Ti{H) is the period of motion in the phase space for the particle 
in the i-th well with energy H: Ti{H) = — ■j=^ where rrii is the width of the i-th well. 

{H{t),K{t)) at each 

. Here i = i{l) is the 



The branching probabilities for the limiting motion Y{t) 



interior vertex O/ are given by p\ 



(0 



.(0 



JO ^ Si) 



-1 



+ cX 



and P2 



(0 



JO 



Si) ^ Si) 



-1 



+ cX 



number of the edge that is above O/ ; p\ \ P2 ^''^ probabilities that the particle enters 
the left edge or right edge that is below O/. The motion inside the edges is described as 
in Corollary (3.1). Since each time when the particle hits the boundary, the perturbation 
of the dynamics is given by independent random variables, the branching at each interior 
vertex is also independent of each other. Therefore we finally have the following theorem: 



>e,5 



(H{X^' ),K{X^' )) on r converges weakly in 



Theorem 3.2 The process Y^ 
uniform topology in the space of continuous functions [0, T] — t- F as first e 4, and 
then 5 ], to the continuous stochastic process Yt = {H{t),K{t)) on F. The process Yt, 
Yq = y(go,Po); inside the edges is described by (3.16). When the process Yt along an edge 
li attached to an interior vertex Oi and situated above Oi reaches Oi, it instantaneously 
leaves Oi and enters edges of F attached to 0[ and situated below Oi on the left and on 



the right with probabilities p\ 



Si) 
'1 



ii) ^ JO ' 



'1 



+ C, 



(0 
P2 



Si) 



Jf+^ 



, respectively. The branching 
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Fig. 9: The case of general potential 



at each interior vertex Oi is independent of each other. 



.(») 



Remark: Let us suppose that 5 = 1: the coefficients of restitution become 1 



(i) 



ec\ — eCu and 1 



ec. 



(i) 






e,l 



respectively. Then the process Y^ ' converges weakly 



(i) (i) 

as e 4, to the process 1^ on F described above with the replacement of c\ and C2 by 



^^) 



Si) 



+ E^(*) and by 



^i) 



.(') 



+ Er/'*\ respectively. 



4 The case of general potential 



The problem mentioned in the beginning of the introduction is considered in this 
section. Let the motion of a particle inside [01,02] be governed by the equation qf = 
—F'{qf), Qq = qq, Qq = pq. When the particle hits the walls, it is reflected with a loss of 

energy. Let H{p, q) = h F{q) be the energy of the particle, p = q. The coefficients 

of restitution at Oj are equal to 1 — eci{H), < e << 1, i G {1)2}, Ci{H) are positive 
and smooth. As was explained in the introduction, we can restrict ourselves to the case 
when F{q) has just one maximum inside [01,02], say, at oq G (01,02) (see Fig.9). Let F 
be the graph corresponding to the Hamiltonian H{p, q) , V~\ = {{p, q) £ M.'^ : ai < q < 02} , 
and y : n — ;• F be the projection of n on F: Y{p, q) = {H{p, q), K{p, q)). 

Denote by X^ = {qf,Pt) nearly-elastic trajectory starting at {qo,Po)- This motion 
has a fast component, which is actually the motion along the elastic trajectories, and the 



slow component Yf 



y(x 



Let Yf 



y/ 



J , and let Y^^, be the continuous piecewise 
linear approximation of Y^, as we considered earlier. As in the model problem, one 
can see that liml^'^ (and liml^^) does not exist. Since we have just two wells, one can 

£^0 £4-0 
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use perturbations of the initial conditions to regularize the problem. Let S,s be a two- 
dimensional random variable with a continuous density which is positive for |x| < 5 and 
is equal to zero if |x| > 5. Denote by Y^' the continuous modification introduced above 
when the initial condition xq = {qo,Po) is replaced by xq+^S: ^t' i^ ^ stochastic process 
since Yq' is random. 

Let us introduce now a stochastic process Yt = {Ht,Kt) on F. First, let Tk{H), 
k £ {1,2,3}, be the period of oscillations of the elastic motion q^ with the initial 
conditions in Y^^{H,K). It is easy to check that 

'^ ' Ja, ^2{H-F{q)) 
n{H) = 2 ' ^ 



T2{H) = 2 



y/2{H - F{q)) 
"2 dq 



_(^H) ^2{H - F{q)) 
Here a±{z) are roots of the equation F{a±[z)) = z, a^{z) < a-^.{z). 
Let Yt = {Ht, Kt) be the stochastic process on F such that Hq = H{qQ,pQ) > H{0), 
Kq = 3; Ht is a deterministic motion inside each edge: 



Ht = -2 



i\Ht){Ht - F{ai'^)) + cf{Ht){Ht - F{af)) 



Tk{Ht 



inside Ik, k G {1,2,3}, where c): {H) = Cj{H), a):' = aj, for /c = 3 and j = 1,2; 
cj^ (H) = Ck{H), aj! = Ofc for k = 1,2 and j = 1, 2. The trajectory Yt hits O in a finite 
rli{qo,Po) Ts(z)dz 

time tn = I — — -— — - — -— — - — -— - since T^iH) ~ ln(H — H(0)) 

Jh(o) 2[ci{z){z-F{ai)) + C2{z){z-F{a2))] '^ ' ^ ^ " 

as H ^ H{0). After hitting vertex O, Yt leaves O immediately and goes to /i or I2 

with probabilities pi = — ,^^,^,, ,^^,„,- and p2 = — ,^^,„,- ttttt^ttt- These 

ci{H{0)) + C2{H{0)) "^ ci{H{0))+C2{H{0)) 

conditions define 1^ in a unique way. 

Theorem 4.1 Under mentioned above conditions, process Y^' converge weakly in 
the space Cqt of continuous functions [0,T] — )■ F, T < 00, as first e J, and then (5 4, 
to the process Yt. 

The proof of this Theorem is similar to the proof of Theorem 2.1 and we omit it. 

One can regularize the problem by introducing random perturbations not of the 
initial conditions, but of the dynamics, as we did in Section 3. Similar to Theorem 3.1, 
one can prove that the process Yt introduced above again serves as the limiting slow 
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motion. These results allow to say that stochasticity is an intrinsic property of the long- 
time behavior of nearly-elastic systems: the limiting slow motion is the same stochastic 
process for different regularizations. 



5 A two dimensional model problem 

Consider a standard billiard model in which a point mass moves freely inside a 
convex, bounded, and simply connected region G in the plane with smooth boundary 
dG (Fig. 10). The orbits of such motion consist of straight line segments inside G joined 
at boundary points according to the rule that the angle of incidence equals the angle 
of reflection. Speed is a constant of motion. The phase space /\ of this system is 
conveniently described as the set of all tangent vectors of fixed length (say, unit length) 
supported at points of the interior of G together with vectors at boundary points pointing 
inward. We parameterize the phase space /\ by the coordinate {g, if), where g G GUdG 
is the position of the particle within G U dG and 99 G [0, 27r) is the angle between the 
positive y-axis and the velocity vector of the particle. The phase space of the motion has 
a global section 11 which consists of all inward-pointing unit vectors on dG. Topologically 
n is a cylinder parameterized by the cyclic length parameter s £ (0, L] along dG (here 
L is the length of dG) and the angle £ [0, vr] between the velocity and the positive 
tangential direction. The Poincare map / : 11 — t- 11 is usually called the billiard map and 
can be described by f{s, 9) = {S{s, 9), Q{s, 9)). Notice that when dG is smooth / is also 
smooth. For a further reference about the billiard model we refer to [9], page 339. 

A classical result in such a billiard model is that / preserves the volume element 
sin 9dsd9 on H. This volume element sin 9dsd9 is usually called the Liouville measure; 
m{s, 9) = sin9 is the density of the Liouville measure. 

We now assume that the particle loses a small amount of energy each time when it 
has a collision with dG. Let c{s, 9) be a smooth function on H, bounded from above and 
below by some fixed positive constants. We assume for brevity that the energy H of the 
particle after the collision with dG at a point x £ H diminishes to maxji:/^ — ec(x),0} 
(rather than to H{1 — ec(x)) as before); here e is a small positive parameter. The 
direction of the velocity vector remains the same as for the elastic system (see Fig. 10). 

We assume that there is another wall of height H{0) separating the region G (see 

const 
Fig. 11). In a finite time of order , e J, 0, each trajectory of the nearly-elastic 

motion, which had at time zero energy Hq > H(0), enters one of the wells 1 or 2 shown 

in Fig. 11 and continues the motion there (we assume that the loss of energy on the 

wall separating the region G is defined in a similar way). Which of the wells is entered, 

in general, depends on e in a very sensitive way: as e 4, the trajectory alternatively 

enters well 1 or 2. As in the previous sections, the nearly-elastic motion has a fast 
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Fig. 10: The billiard model 

and a slow components. In the case of one degree of freedom, the fast motion has just 
one normalized invariant measure, and a unique, independent of the initial conditions, 
limiting slow motion inside the edges exists without any regularization. This is not 
the case now. Besides the Liouville measure, other invariant measures on the Poincare 
section 11 exist; for instance, there are periodic points on dG. On the other hand, it is 
natural to assume that the reflections on dG are undergoing small random perturbations 
of intensity 5 « 1 so that the resulting motion is a stochastic process depending on 
two small parameters e and 5. We choose a class of random perturbations in such a 
way that the perturbed system induces on 11 a unique invariant measure, namely, the 
Liouville measure (compare with [5]). It turns out that the perturbations in the fast 
motion which we introduce provide a regularization of the slow component not just 
inside of the edges of the corresponding graph but also near the vertices, so that the 
weak limit of the slow component of the regularized motion (after an appropriate time 
change) exists on the whole graph as first e and then 5 tends to zero. Probably, the 
class of permissible random perturbations leading to the same limiting slow motion can 
be essentially extended. But one should keep in mind that the existence of a rich set of 
invariant measures for the billiard map shows that the limiting slow motion calculated 
in this section has, in a sense, a restricted universality: other regularizations of the fast 
motion may lead to different slow motions. 

We regularize this problem by considering the following scheme of small stochastic 
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perturbation of the system. Consider an one-dimensional diffusion process if on [0, vr] 
starting from point 9. The diffusion process is governed by a second order ehiptic 
operator C which is self-adjoint with respect to the Liouville measure. Such an operator 

can be written as Cu = 7:^(a(0)—r;:), where a is a smooth positive function of 6. 

We assume that the diffusion process if has instantaneous reflection at the boundary 
points ^ = 0, TT. Notice that there is a singularity of our operator C at the boundary 
points ^ = 0, vr. Inspite of this singularity the boundary points ^ = 0, vr are accessible 
and therefore we need boundary conditions: we choose instantaneous reflection at these 
points. 

Choose a small constant 5 > 0. Let the position of the particle after one collision 
be X = (s,0) G n with an energy H > 0. Then the next position on 11 will be not 
exactly x but Z^ = (s, /^), and the energy will diminish to H — ec{x) (see Fig. 11). And 
the subsequent motion starts from this new point Z^. The same kind of perturbation is 
repeated again and again for each collision, independent of other collisions. This scheme 
of perturbation models the fact that each time when a collision happens, the particle 
changes its direction of motion a little bit, due to the random perturbations. 

After adding the small stochastic perturbation, the evolution of energy and the mo- 
tion of the particle is described by a continuous time Markov process N^' = {H^' , A^' ) 
on [0, +oo) X /\. The energy Hp changes just on the boundary: if A^' = x = {s,6) £ H, 
the energy instantaneously decreases on ec{x); Hi' is right continuous. The process 
A^' = ((7j ' , ifl' ) has a piecewise linear first component: it moves uniformly along the 
chords connecting successive points where gl' hits dG with the constant speed v/2i7^' . 
Let M^' = [Hi' ,gl' ) (a sample trajectory of M^' is the thin line in Fig. 11). The tra- 
jectory of M^' is right continuous and has limit from the left. We define a continuous 
modification M^' = (Hi' , g^' ) of M^' by taking H^' as a piecewise linear modification 
of H^' (a sample trajectory of M^' is the bold line in Fig. 11). Let N^' = {Hi' , A^' ). 

Consider also a jumping Markov chain {X„}„>i on IT, defined as follows: from a 
point X^ = (s„,^„) G n the chain goes in one time unit to Xf^^^ = (Sn+i, ^n+i) = 

Suppose the particle of unit mass starts its motion from a point x = (s, ^) G IT and 
the energy of the particle is H. Let T(x, H) be the time that the particle starting from 

x first reaches dG. Let l{x) be the chord starting from x G IT (see Fig. 10). Let L{x) be 

Lix^ 
the length of l{x). We have T(x, H) - 



V2H 
We identify points of each well having the same energy. The set obtained after 

such an identification is a graph F; the metric on F is defined by the distance along the 

edges of F. This graph has one interior vertex O corresponding to the energy level H{0) 

connected with three edges: I3 corresponding to trajectories with energy greater than 

H{0), /i corresponding to trajectories with energy less than H{0) situated within well 
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x=f(xo)=(s,e) 



H-sc(x) 



H(0) 




Fig. 11: The billiard model with energy loss and perturbation 



1 (to the left of O), and I2 corresponding to trajectories with energy less than H{0) 
situated within well 2 (to the right of O). 

We consider the process Y^' = {Hi' , K^' ) on the tree F, where K^' is defined 
as follows: K^' = 3 if the energy of the particle is greater than H{0), K^' = 1 if the 
particle has energy less than H{0), while A'^''' is in the potential well 1, and K^' = 2 
if the particle has energy less than H(0), while N^' is in the potential well 2. The 
process Y^' is the slow component of N^' as e J, 0. We show that the rescaled process 
Yf^,' converges, as first e 4, then (5 4, 0, to a stochastic process Yj on the tree T. Such 
a stochastic process has deterministic motion within each edge of the tree F and has 
stochasticity only at the interior vertex O. 



Let us first state some auxiliary results: 



Lemma 5.1 The diffusion process if on [0,7r] governed by the operator Cu - 

;:-r;7(a(^)^) with instantaneous reflection at boundary points 9 = 0,tt, has — sin( 

2sin6'd6'^ ^ ^ d9^ ■' ^ ^ ' ' 2 

as its unique invariant density. 



Suppose small r > 0, let n(r) = {x = {s,6) £ Il;r < 6 < n — r}. 



Lemma 5.2 For any 6 > 0, for any small enough r > 0, the Markov chain 
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{X^}n>i on n(r) satisfies the Doeblin condition, and has only one ergodic compo- 
nent on n(r). The invariant measure for the process {X„}„>i on H has a density 

d(x) = — m(x) = — sinO. 
^ ^ 2L ^ ' 2L 

Lemma 5.3 Within edge I3 of the graph T, linilinii^w = Ht in probability, where 
Ht satisfies the following differential equation: 

c{x)m{x)dx 

, Ho = H'/ (5.1) 



dHt _ JJn rr _ TTe,S 

-df-~ ~ ,Ho-H, 



T{x,Ht)m{x)dx 
n 

In exactly the same way as this lemma one can have similar results for limiting 
deterministic motion within edge /i and 12- We omit details here. 



Lemma 5.4 Denote the area of the region hounded by the convex curve in Fig. 10 
by A. Then 

L{x)m{x)dx = 2ttA 



/n 

where m{x) is the density of the Liouville measure on U and L{x) is the length of the 
chord l{x). 

Using Lemma 5.4, equation (5.1) can be simplified as 
And (5.2) has an explicit solution 



H,= (^„-^-ff .(.)„.(:■:)..) \0<t< ^'-^'f^ - V^^' . (5,3) 
V 2V27,AJJt, } n ^(^)„(3.);j^ 

Lemma 5.5 For any e > 0, any 5 > the invariant measure for the process A^' 
is proportional to the Lebesgue measure on /\. 

Lemma 5.6 For any 6 > 0, 

limPjA'^''' finally falls into well 1} = limPjA^^''' finally falls into well 1}, 
sio £4,0 

limPjA'^j''' finally falls into well 2} = limPjA^^''' finally falls into well 2}. 
e4,o £4,0 
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At the beginning of this section we have already defined the process y/' = {Hi' ,Kl' ) 
on the tree T. Let IIj, i G {f , 2} be the set of those points on 11 which correspond to 
well i , i £ {1,2}, respectively (see Fig. 11, but notice that it only shows the s coordinate 
in the space H). 



Lemma 5.7 



limlimPjA'"^^' finally falls into well 1 } = 

54-0 e4-0 



c{s, 9)m{s, 6)dsd9 



c{s,9)m{s, 9)dsd9 
n 

c{s, 9)m{s,9)dsd9 



limlimPjA^j'^' finally falls into well 2 } = 
Sio £4,0 



c{s, 9)m{s, 9)dsd9 
n 



Now similarly as in Section 2 and Section 3, we consider a piecewise linear modifica- 
tion of the process Y^' , and denote such a process by Y^' . We have Y^' = {Hi' , Kl' ). 
For all e > 0, for all < t < T, for some uniform constant C > 0, P{\Y^' — Y^' \ < 
Ce) = 1. We define a stochastic process Yj on F as follows: it has deterministic motion 
within each edge of the tree F and only has stochasticity at the interior vertex O. Let 
Yt = {H{t),K{t)) where H{t) is the energy of the particle and K{t) is the number of the 
edge of the graph F. The deterministic motion within each edge of the tree F is defined 
as follows: on edge 3 the first component H{t) of Yt satisfies 

Ht = f /ffo /t— ff c{x)m{x)dx\ , < t < to , (5.4) 

and on edge 1 and 2 the differential equation becomes 

Ht = LjH{0) - -^^^ Ij c{x)m{x)dx] , to < t < t, , (5.5) 

Ht = (VH{0) - ^~^° // c{x)m{x)dx\ , to<t<t2, (5.6) 

where Hi and 112 are defined as before, and Ai, A2 are the areas of the domains in G cor- 
responding to well 1 and 2. The process Yj, starting from a point Yq = {Hq, 3), will reach 

, . . ^. . 2TTA{./m'Q-J2H{0)) , .„. 

the interior vertex (J m time to = F7 < 00 and will instantaneously 

/ / c{x)m{x)dx 

I I c{x)m{x)dx 

leave O and enter one of the edges 1 or 2 with probabilities pi = — j-jr 

/ / c{x)m{x)dx 
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c{x)m(x)dx 
and p2 = " I- ^ respectively. After Yt enters edge 1 or 2 it will move 

c(x)vn(x)dx 
'n 
deterministically according to the given differential equation (5.5) or (5.6) till time 

ii = io H — FT ) ^ £ jli 2|, respectively, when it hits energy level zero. 



c{x)m{x)dx 
Hi 

These lemmas imply the following: 



>e,5 



Theorem 5.1 The process Y / converges weakly, as first e | then 610, to Yt 



Let us prove now these lemmas. 

Proof of Lemma 5.1: Let r be the first time for process if to exit from the interval 
[0, vr]. By considering u{6) = E^r as the solution of the corresponding Dirichlet problem, 
one can check that u{9) is finite so that the boundary points = 0, vr are accessible. 

Now we consider another diffusion process if which starts from point 6 and is gov- 

erned by the operator Cf = ~~t^{(^{^)~t^) with instantaneous refiection at the boundary 
points = 0, TT. The invariant denisty for the process if is the uniform distribution with 
the density — . 

Our process if can be obtained from if by taking a time change dt = sin Odt. There- 
fore invariant measure for if is — sin^d^ (see, [8], Chapter 5 for a reference concerning 
the random time change). ■ 

Proof of Lemma 5.2: For any small r > 0, any 5 > 0, any ^i, ^2 G V"-, t^ — f]^ we have 
Pi{5; 01,62) > a{r) > for a positive constant a{r) > depending on r, where pi{-; •, •) 
is the transition density for the process if. 

Let r > 0. Let n(r) = {x = (s,^) G n;r < ^ < vr — r}. Since the function 
Q{x) = Q{s,6) is continuous for x £ n(r) and the set n(r) is compact, Q(x) must attain 
its maximum and minimum on n(r). It is easy to see that there exist A(r) > such that 
min(r, vr — max Q(x), min Q(x)) = A(r) > (otherwise the particle must start from 

the tangential direction and keep on moving along the boundary, which contradicts the 
fact that 6 £[r,TT - r]). This implies that /(n(r)) C n(A(r)) with < A(r) < r. 

We check the following form of Doeblin condition: For any x G n(r), there exist 
an integer uq > 1 and a positive e > such that for any Borel measurable function A 
in B{Il{r)), iy{A) < e, we have p^"''''{x,A) < 1 — e. Here i^(-) is the standard Lebesgue 
measure on H. 

Let x = (si,^i) G n(r). Let (53,^2) = (5(si,^i), e(si,^i)) G n(A(r)). Let a 
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point (53,^3) G n(r), S3 / 82- Let C("525-53) G (0, vr) be the unique angle such that 
Sis2,Cis2,S3)) = S3. Let the set W = {S{s2,0);e £ [A(r),7r - A(r)]}. We choose 
£1 > such that iy{U{r)) - ei > £2 > 0. Then for iy{A) < ei we have iy{U{r)\A) > 
£2- We first consider the case when iy{{W x [r,9 — r]) n {Il{r)\A)) > £3 > 0. This 
imphes that there exist £4 > such that nies({s3 G W; (53, ^3) G Il{r)\A for some ^3 G 

[r,TT — r]}) > £4 > 0. Here nies(-) is the standard Lebesgue measure on M. We have, 

dC 

> £5 > for S3 G W. 



for some £5 = £5(A(r)) > depending on A(r), that 

dS3 

Let K = {Cis2, S3); S3 G M^, (53,6*3) G Il{r)\A for some ^3 G [r,TT - r]}. We see that 
mes(-ftr) > £4£5 > 0. We also notice that K C [A(r),7r — A(r)]. Let us define the set 
M = {6*3 G [r, vr - r]; (S3, ^3) G Il{r)\A for some S3 G W}. There exist £5 > such that 
mes(M) > £6 > 0. Now we have 

p(-^\x,U{r)\A) 

Pi{6;e2,e)pi{d;e{s2,e),e3)ded93 



leaGM JeeK 

> [ [ a(A(r))a(A(A(r)))d0d^3 

> a(A(r))a(A(A(r)))£4£5e6 > . 

Therefore p^'^\x,A) < 1 — a(A(r))a(A(A(r)))£4£5£6. So in this case we can choose 
no = 2 and < £ < min(a(A(r))a(A(A(r)))£4£5£6,£i). 

Now let us consider the case when ^{{W x [r,9 — r]) n (n(r)\^)) = 0. Since we 
have u{U{r)\A) > £2 > 0, we have iy{{{[0, L]\W) x [r,e - r]) n {U{r)\A)) > £2 > 0. 
In this case we can take one more transition from (s2,^2) to some other set of positive 
measure lying away from [0, L]\H^ (provided that r is small enough) and carry out a 
similar argument as above. We conclude in this case that no = 3. 

It is easy to see that this chain has only one ergodic component n(r), since from 
any point x G n(r) the chain can come to any neighborhood of any other point y G n(r) 
in several steps. 

Since m{x) is an invariant density for both the billiard map / and the diffusion 
process if, it is also invariant for the chain {X^}n>i. Therefore, after normalization, 

the unique invariant density for the chain {X„}„>i is d{x) = —-m{x). ■ 

~ Zlj 

In particular. Law of Large Numbers holds for this Markov chain {X^}n>i on 11. 

Proof of Lemma 5.3: Suppose there are n = n{e) collisions happened during 

time [-, ), and the corresponding successive positions of the particle on 11 are 

X^j^^, ...,X^^^. It is clear that n — ;■ 00 as £ | and ne ~ 0(At) as £ | 0. For small 
At > 0, we have 



30 



rre,S rre,S .^^ At ^ c{X',+,) + c(X^+,) + ... + c(X,\J 



We also have 

At _ tinitial - t final + T{Xf^_^_^, Hi) + T{Xi^^^, H2) + ... + T{Xj^^^^, Hn) 

ne n 

where Hi , 1 < i < n are the energies of the particle at the i-th collision (before it loses 
some energy); tiniual is the time between the starting time - of the motion and the time 

of the first colhsion with dG; tf,r.ai is the time between the ending time ^ of the 
motion and the time after that when the particle has next (the n + 1-th) collision with 
the boundary of G. Since function T{x, H) is Lipshitz in H (when H > h > 0) with 
respect to x, we have \T{x, Hi) - T{x, H^^Jl < G\Hi - H^^J ~ 0(At) as e ^ 0. Therefore 
from (5.2) , by using the Law of Large Numbers for this Markov chain , we get 



At 
ne 



ET{X,Ht/^) + 0{At) 



in probability, as e J, 0. 

By using the Law of Large Number for our Markov chain we also see that 

c{xi^i)+c{xi^,) + ...+c{xi^j ^ 

n 

in probability, as n — t- oo. 

Here the random variable X has the stationary distribution of this chain on H with 
density d{x) = —-m{x). Therefore 

^^H(t+At)/e-Ht/e _ MX) 



40 At ET{X,H^^J + 0{At) 

in probability. 

This means that as e J, 0, the rescaled process H^', within edge Is of the tree T 
converges in probability to the trajectory of the deterministic process defined by the 
differential equation 



^fj6 II c{x)m{x)dx 



T{x, Ht)m{x)dx 
in 



Notice that the this is independent of 6. Therefore letting 6 ^0, we get the desired 
result. ■ 
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Fig. 12: A problem in integral geometry 

Proof of Lemma 5.4-' This is a known result in integral geometry (see [11]). But 
we will still provide a short proof here in order for the reader to understand some 
calculations later in this section. 

Let the equation of the line intersecting the closed curve dG be written in the form 
(see Fig. 12) 



Let the curve dG be 



X cos if + Y sin (p = p . 



X = X{s),Y = Y{s) 



where s is the arc length parameter. 

Let r be the angle of the tangential direction at point {X, Y), we have 



vr 



COST = X'{s) , sinr = Y'{s), — + ip — t 



We then have 



dip = dO + dr = dO + nds . 
Here k = k{X,Y) is the curvature at point {X,Y), and 
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dp = X'{s) cos ifids — X{s) sin ipdip + Y' [s) sin (pds + Y{s) cos ipdip 

= {X'{s) cos if + Y'{s) sin (/3)(is — X{s) sin ip{nds + d^) + l^(s) cos {p{nds + d^) 

= (X'(s) cos (p + y (s) sin (y9 + n(—X{s) sin (/? + Y{s) cos ip))ds + (— X(s) sin (/? + l^(s) cos ip)d9 

These equalities imply that 

dp A dip 

= {X' {s) cos 99 + Y'{s) svnip + k{—X{s) sine/? + l^(s) cos ip))ds A d^ 

+k(— X(s) sin if + y (s) cos 99)^^ A ds 
= {X'{s) cos tp + y (s) sin (/9)ds A dO 
= (cos r cos ip + sin r sin ip)ds A d^ 
= cos(r — ip)ds f\d9 = svnOds A dO . 

So we have that 

/ / L{s,e)smedsAde = 2 / L{s,e)dp A dip = 2ttA. 
Jo Jo Jo Jo 

We have used the fact that / L{s, 9)dp = A, where P is the height of G projected 

Jo 
onto the direction perpendicular to l{s,9) (recall that l{s,9) is the chord starting from 

{s,e)eu). u 

We would like to mention that the relation dp A dip = sin Ods A dO will be used in 
some of our later calculations. In future calculations, we will use the abbreviated form 
dpdip = sin 9dsd9. 

Proof of Lemma 5.5: From Lemma 5.2 we already know that the embedded chain 
{X^}n>i has a unique invariant measure on 11. Since our process A^' is purely de- 
terministic within /\ \n, and any point in /\ \n can be reached in finite time from a 
point on H, A^' must be positive recurrent on /\ and therefore it has a unique invariant 
measure on /\. 

We use the Khasminskii formula (see, [6], on page 185, formula (4.1)) for the ex- 
pression (up to a constant factor) of the invariant measure /!(•) for A^' on /\ through the 

invariant measure —--m{x)dx for the embedded chain {X„}„>i on H: for any measurable 
set B C /\, 

KB) = ^jf m{x)dxW.^ j ' XB{K'^)dt . (5.8) 

Here rf is the first time when A^' , Aq' = x G II, reaches 11 again. Note that for a 
given X, rf as well as the motion A^' , < t < rf , and in particular, the time spent by 
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A^' in B, are not random, so that the expectation sign in (5.8) can be omitted. Since 
the invariant measure of our imbedded chain and of the bihiard map are the same, we 
conclude from (5.8) that Jl{B) coincides with the invariant measure of the non-perturbed 
bilhard. As is known, the Lebesgue measure is invariant for the bihiard. Thus fl{B) is 
the uniform distribution on /\. ■ 

Proof of Lemma 5. 6: This follows from the fact that 

P{ max \Hf^ - Hf'^\ < Ce} = 1 

0<t<T<oo 

where C > is a constant. ■ 

Proof of Lemma 5. 7: Consider a time change in the process A^^ ' . If at time t, A^' 

~ Lix) 
lies on a chord lix) starting from x = is, 9) G IT, then we set dt = dt. Here 

c{f{x)) 
f{x) = f{s, 9) = {S{s, 9), @{s, 9)) is the billiard map and L{x) is the length of the chord 

l{x). The density of the invariant measure for A~ (with respect to Lebesgue measure) 

cifix)) 
is proportional to — . Here x = x{g, (p) is a point in H such that (^f, 99) G A lies on 

the chord l{x) starting from x. This time change corresponds to multiplication of the 

generator of the process Nf' by , „ , . , . 

4f{x)) 
The iZ-component of the process obtained after time change moves down uniformly 

with a non-random speed, so that it hits level H{0) at a non-random time t = t^Ho). 

Let Ai = {(5, 99) G A : f{x(.9,^)) ^ Hi}, A2 = {{g,ip) G A : f(.x{g,ip)) G U2}. Here 

ni,n2 are defined as before, i.e., Hj, i G {1,2} is the set of those points on H which 

correspond to well i , i G {1,2}. 

By using Lemma 5.5, we see that 

lim P{N~ finally falls into well 1 } 
JJJai L{x{g,ip)) 

= D ff ^H^P^ \ r ^'^ dLis, 9) 1 sin 9dsde 

JJf'H^i) ^(s>6') Wo / 



Dj c{f{s,9))m{s,9)dsd9 



= D c{s,9)m{s,9)dsd9. 

Here m{s,9) = sm9 is the density of the Liouville measure. We have used the fact 
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that dgdip = dL{s, 6)dpd(p = dL{s, 9) sinOdsdO (notice here that dg = dLdp), as was cal- 
culated in Lemma 5.4, and the invariance of Liouvuille measure under the billiard map. 

Normalizing constant D = I fff -^j^^j^^^^dgd^ j = ( // c(s, e)m{s, 0)dsdej . 

11 c{s,e)m{s,e)dsde 
Similarly, we have limP{A^~' finally falls into well 2 } 

^■^° * " c{s,9)m{s,9)dsde 



^ _ 'n 

The trajectory of N~ is the same as N^' . Taking into account this fact and 

Lemma 5.6, the statement follows. ■ 
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